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ABSTRACT 



■ The temperature (T) and entropy (S) fields of baryonic gas, or intergalactic 

O ■ medium (IGM), in tlie ACDM cosmology are analyzed using simulation samples 

^ ■ produced by a hybrid cosmological hydro dynamic /A^-body code, effective cap- 

turing shocks and complex structures with high precision. We show that in the 
nonlinear regime the dynamical similarity between the IGM and dark matter will 
be broken in the presence of strong shocks in the IGM. The heating and entropy 
production by the shocks breaks the IGM into multiple phases. There are no 
^ ■ single-value relations between T or S and the mass densities of the IGM (pigm) 

^ ■ or dark matter (pdm)- The probability distribution functions of the temperature 

and entropy fields are long-tailed, with a power law decreasing on the sides of high 
temperature or high entropy. These fields are therefore intermittent. The mean 



CD 



in 
o . 

'nI" . entropy, or the cosmological entropy floor, is found to be more than 100 h 

keV cm^ in all regions when z < 1. At redshift z ~ 2 — 3, high-entropy gas 
P_|. (5* > 50 h~^^^ keV cm^) mostly resides in areas on scales larger than 1 Mpc 

Q . and with density pdm > 10^ (in units of 'pdm)- Therefore, gravitational shocks are 

an effective preheating mechanism of the IGM and probably enough to provide 
. the entropy excess of clusters and groups if the epoch of the gas falling in cluster 

cores is not earlier than 2; ^ 2 — 3. On the other hand, at redshifts z < 4, there is 
^ . always a more than 90% volume of the low dark matter mass density {pdm < 2) 

regions filled by the IGM with temperature less than 10^'^ K. Therefore, the 
multiphased character and non-Gaussianity of the IGM field would expain the 
high-temperature and high-entropy gas observed in groups and clusters with low- 
temperature IGM observed by Lya forest lines and the intermittency observed 
by the spikes of quasi-stellar object's absorption spectrum. 
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1. Introduction 

The standard scenario that provides the theoretical skeleton for the formation of cosmic 
large-scale structures in the universe is the hierarchical clustering scenario. This scenario 
assumes that massive dark matter halos are formed via the processes of gravitational col- 
lapse, merging of smaller structures, and accretion. Light-emitting and -absorbing objects 
form when the baryonic gas, i.e., the intergalactic medium (IGM), falls, cools, and condenses 
around and within the dark halos. The evolution of the IGM in this model is largely deter- 
mined by the underlying dark matter. It is usually assumed that the dynamical behavior of 
the IGM field can be obtained from the dark matter field by a similarity mapping between 
the two (e.g.. Kaiser 1986). The perturbations in the density and velocity distribution of 
the IGM are considered to be the same as those of the dark matter field point by point on 
scales larger than the Jeans length scales. 

In the linear regime, this similarity mapping is correct. Even when the perturbations 
of the initial IGM distribution are different from those of the dark matter, linear growth 
modes will lead to similarity on scales larger than the Jeans length (Bi et al. 1992; Fang 
et al. 1993; Nusser 2000; Nusser & Haehnelt 1999). However, in the nonlinear regime, the 
similarity between the IGM and dark matter fields on scales larger than the Jeans length 
may not always be present. It has been shown in hydrodynamics that a passive component 
generally decouples from the underlying field during nonlinear evolution (for a review, see 
Shraiman & Siggia 2000). For instance, a passive component might be highly non-Gaussian, 
even when the underlying field is Gaussian (Kraichnan 1994). This nonlinear decoupling is 
common in stochastic systems consisting of a "passive component" in an underlying mass 
field. 

In the system consisting of the IGM and dark matter, the IGM can be considered a 
"passive component", because the evolution of the IGM is driven by the gravity of the 
underlying dark matter. Moreover, the system of the IGM and dark matter is stochastic 
because of the randomness of the initial perturbations of the cosmic mass and velocity 
fields. The mass and velocity fields of the IGM and dark matter are random variables. The 
importance of the stochastic nature in the evolution of cosmic mass and velocity fields has 
been emphasized in many works (e.g., Berera & Fang 1994; Jones 1999; Buchert et al. 1999; 
Matarrese & Mohayee 2002; Coles & Spencer 2003; Ma & Bertschinger 2004). Therefore, in 
the nonlinear regime, we expect that the similarity mapping between the IGM and the dark 
matter may no longer work. 

The dynamical mechanism of violating the similarity between the IGM and dark matter 
has been discussed in the early study of structure formation (Shandarin & Zel'dovich 1989). 
As the dark matter particles are coUisionless, the velocity of the dark matter particles will be 
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multi valued at the intersection of the dark matter particle trajectories. On the other hand 
the IGM, being an ideal fluid, has a single-value velocity field. Thus, at the intersection of 
the dark matter particle trajectories, the density of gas would be discontinuous, and shocks 
or complex structures occur. The development of shocks or complex structures breaks the 
similarity between the IGM and dark matter and thus leads to the statistical discrepancy 
between the IGM distribution and dark matter mass field. 

Strong shocks and complex structures result in the heating and the entropy increase 
of the IGM. The cosmic baryonic gas involved in strong shocks should hence undergo an 
entropy increase. The dark matter remains unaffected. Consequently, the IGM will be 
in thermodynamically multiple phases. The relations between the temperature and mass 
density, or entropy and mass density of the IGM will not be described by an single-value 
relation and will vary from place to place. It would be interesting to study the multiphased 
properties of the IGM caused by the statistical discrepancy between the fields of the IGM 
and dark matter. In this paper, we focus on the temperature and entropy fields of the IGM. 

Aside from theoretical interests, this study is also motivated by the problem of the so- 
called entropy excess, i.e., a high entropy fioor is observed in nearby groups and low-mass 
clusters (Ponman et al. 1999; Lloyd-Davies et al. 2000). Entropy excess indicates that 
some processes that violate the dynamical self-similar scaling are important in structure 
formation. To explain this result, many models for IGM preheating have been proposed, 
such as photoionization heating (e.g., Efstathiou 1992; Gncdin 2000), supernova (SN) explo- 
sions, active galactic nucleus (AGN) activities (e.g.. Cole et al. 1994; Somcrville & Primack 
1999), gravitational shocks associated with clusters (Cavaliere et al. 1998) or warm- hot gas 
(Valageas et al. 2003). 

The decoupling of the IGM and dark matter in the nonlinear regime inevitably leads to 
a dynamical break of the similarity between dark matter and gas. Hence, the temperature 
and entropy fields of the IGM formed by this evolution should be treated as the background 
or cosmological "floor" . The entropies produced by other mechanisms are superpositions on 
this background. Many previous studies assumed an uniform heating of the IGM all over 
the space (e.g.. Mo & Mao 2002). However, this assumption will not hold if the background 
temperature and entropy are inhomogeneous and multiphased. Some studies of the grav- 
itational shock heating have assumed that the heating and entropy increase occur only in 
nearby groups and clusters. This is to assume that shock heating occurs mainly in high- 
density regions. However, the decoupling evolution and gravitational shocks can happen in 
high mass density regions (hosting groups and clusters), as well as in low mass density areas. 
Thus, shocks can also heat the IGM at low mass density areas. 

The last motivation for this study is the progress made in the art of numerical algorithms 
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for hydrodynamics over the last few years. The weighted essentially nonoscillatory (WENO) 
scheme (Fedkiw et al. 2003; Shu 2003) has proved to be significantly superior over the 
piecewise smooth solution algorithm containing discontinuities. It has been successful in 
capturing shocks and complex structures and in calculating the entropy in both shocked 
and unshockcd regions. Therefore, it would be interesting to calculate the temperature and 
entropy fields of the IGM with this code. 

We analyze the temperature and entropy fields of the IGM in the ACDM cosmology with 
simulation samples produced by a hybrid hydrodynamic/A^-body code. This code effectively 
simulates shocks and complex structures with high precision. The paper is arranged as 
follows. The basic equations and the mechanism of the IGM-dark matter discrepancy are 
addressed in §2. In §3 we present the hydrodynamic cosmological simulation scheme. In §4 
we show the basic properties of the temperature field of the IGM. The entropy field of the 
IGM is calculated and analyzed in §5. Finally, in §6 we present our conclusions followed by 
a discussion of the results. 



2. The basic equations 

2.1. Hydrodynamic equations for the IGM 

The IGM is assumed to be an ideal fluid with polytropic index 7. The hydrodynamic 
equations for the IGM in the expanding universe can be written as 

U + d.f{U)=f{t,U) (1) 

where di = d/dX'- (i = 1,2,3); X* denote the proper coordinates, which arc related to 
comoving coordinates by X* = a{t)x^, a{t) being the scale factor. The quantity U in 
equation(l) contains five components: 

U^{p,p,v,E) (2) 

where p is the comoving density of the IGM, v = {vi} is the peculiar velocity, E = P/{'-f — 
1) + |pv^ is the total energy per unit comoving volume, P = a^p, and p is the pressure of 
the IGM. The quantities f\U) in equation(l) are given by the conservation laws of mass, 
momentum, and energy as 

f{U) = [pv\pv'v\p{v-'f + P,pv%\v\E + P)\ 

f{U) = \pv\pv\\pv\\p{v^f + Py{E + P)] (3) 
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The "force" term f{t, U) on the right-hand side of equation(l) is given by 

f{t, U) = (0, --pv + pg, -2-E + pv • g - A,„,). (4) 
a a 

The term of —{a/a)pv is from the expansion of the universe. The term of Arad in equation 
(4) is given by the radiative heating-coohng of the baryonic gas. The gravitational force 
g = — V0 is produced by the coUisionless dark matter, given by 

V^(/) = AnGa'^ptotStot, (5) 

where the operator V acts on the comoving coordinate x; dtot = [Ptot(x, t) — Ptot]/ptot, and 
Ptot is a sum of the comoving baryon and dark matter mass density. The total mean density 
is ptotit) = l/GnGt^ oc a^'^. The gravitational potential 4> is zero (or constant) when the 
density perturbation 6tot = 0. 



2.2. The statistical discrepancy between the IGM and dark matter 

Before embarking on the numerical calculations, we point out some qualitative features 
of the IGM evolution given by equation (1). To sketch the evolution of gravitational clus- 
tering, we ignore the heating and cooling. Consider the case in which all thermal processes 
are described by the polytropic relations p(x, t) oc p"' (x, t), T oc p'^~^, or T = To(l -|- Sy^^, 
where S = [p(x, t) — p]/p is the IGM mass density perturbation. Prom equation (1), the 
momentum equation is 

dav 'ykeT V5 

-— --F(V-V)V = 7" rr - V0 (6) 

dt ^ ' prup (1 + 5) ^ ^ 

where the parameter p is the mean molecular weight of the IGM particles and rUp is the 
proton mass. 

For growth modes in the perturbations, velocity is irrotational. We can then define a 
velocity potential by 

V = --Vif. (7) 
a 

In the linear approximation of the temperature-dependent term, equation (6) yields 



where the coefficient u is given by 



pmp{d\n. D{t) / dt) ' 



(9) 
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in which D{t) describes the hnear growth behavior. The term with u in equation (8) 
acts hke a viscosity (due to thermal diffusion) characterized by the Jeans length kj — 



Since the gravitational potential (f) is mainly given by the random density field of dark 
matter, the term of equation (8) plays the role of a random force. The nonlinear equation 
(8) is exactly the stochastic force-driven Burgers equation or the Kardar-Parisi- Zhang (KPZ) 
equation (Kardar et al. 1986; Berera & Fang 1994; Barabasi & Stanley 1995). The behavior 
of the field ip given by equation (8) depends on the Reynolds number, defined by 7^ = 
(rc/flo)"'^'^ (e.g., Lassig 2000), where Vc is the comoving correlation length of the random field 
(f) and ao is the dissipation length. In the case of equation (8), the dissipation length is given 
by the Jeans smoothing equation (9), or ao = i^^/^ry^(0^)~^/^((ilnD/(it)~^/^. Thus, we have 



where Sdm{k) is the Fourier component of the dark matter density contrast for the wavenum- 
ber k. To derive equation (10), we assume that the gravitational potential (p is only given 
by the dark matter mass perturbation. 

It is well known that the Burgers equation does not reduce, but increase initial chaos 
(Kraichnan 1968). When the Reynolds number (eq. [10]) is larger than 1, the (p field is 
in Burgers turbulence (e.g., Lassig 2000). That is, when the Reynolds number is high, an 
initially random field always results in a collection of shocks with a smooth and simple 
variation of the field between the shocks. 

The correlation length Tc of the gravitational potential is larger than the Jeans length, 
i.e., kjTc > 1. Therefore, TZ can be larger than 1 on scales larger than the Jeans length 
{k < kj), even when 6dm{kj) is of the order of 1. This means that the Burgers turbulence 
on scales larger than the Jeans length can occur in the IGM ip field, when the dark matter 
mass density perturbations are still quasi-linear or weakly nonlinear {Sdm{k) ^ 1) on those 
scales. This results in the decoupling of IGM velocity field from that of the dark matter 
on scales larger than the Jeans length. When the Burgers turbulence develops in the IGM, 
the field (f is characterized by strong intermittency and contains many discontinuities, or 
shocks. The probabihty distribution function (PDF) of cp is long-tailed (Lassig 2000). That 
is, the randomly distributed shocks and intermittent spikes lead to the statistical discrepancy 
between the velocity fields of the IGM and dark matter. 



{ayt^){ump/^kBTo). 




(10) 
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3. Cosmological hydrodynamic simulations 

3.1. Method of hydrodynamical simulations 

Prom the result of §2, it is clear that we need to focus on shocks and discontinuities in 
our simulations. In the system of the IGM gravitationally coupled with dark matter, the 
temperature of the IGM generally is in the range 10^ — lO^K, and the speed of sound in the 

IGM is only a few km to a few tens km s~^. On the other hand, the rms bulk velocity 
of the IGM is of the order of hundreds km s^^ (Zhan & Fang 2002). Hence, the Mach 
number of the gas can be as high as ~100 (Ryu et al 2003). The shocks and discontinuities 
are extremely strong. Hence we will not use a Lagrangian approach such as the smoothed 
particle hydrodynamic (SPH) algorithm. A main challenge to the SPH scheme is the handling 
of shocks or discontinuities, because the nature of SPH is to smooth the fields considered 
(e.g. Borve et al. 2001; Omang et al. 2003.) 

We take a Eulerian approach instead. A well-known problem of the Eulerian algorithm 
are the unphysical oscillations near a discontinuity. An effective method to reduce the 
spurious oscillations is given by designed limiters, such as the Total- Variation Diminishing 
(TVD) method (Harten 1983) and piecewise parabofic method (PPM; GoUella & Woodward 
1984). A problem of the TVD method is that the accuracy necessarily degenerates to first 
order near smooth extrema (Godlewski & Raviart 1996). This problem will cause errors 
in calculating the temperature and entropy changes because they are determined by the 
difference of the thermal energy P/(7 — 1) on the two sides of the shock. When the Mach 
number of gas is high, the thermal energy P/(7 — 1) is very small compared to the kinetic 
energy pv^/2. To overcome this problem, the essentially nonoscillatory (ENO) and WENO 
algorithms were proposed (Harten et al. 1986; Shu 1998, 2003; Fedkiw et al. 2003). They 
can simultaneously provide a high-order precision for both the smooth part of the solution 
and sharp shock transition (Liu et al. 1994; Jiang & Shu 1996). 

The WENO scheme WENO has been successfully applied to hydrodynamic problems 
containing shocks and complex structures, such as shock- vortex interaction (Grasso & Piroz- 
zofi 2000a, 2000b), interacting blast waves (Liang & Ghen 1999; Balsara & Shu 2000), 
Rayleigh-Taylor instability (Shi et al. 2003), and MHD (Jiang & Wu 1999). WENO has also 
been used to study astrophysical hydrodynamics, including stellar atmospheres (Del Zanna 
et al. 1998), high Reynolds number compressible flows with supernova (Zhang et al. 2003), 
and high Mach number astrophysical jets (Carrillo et al. 2003). In the context of cosmo- 
logical applications, WENO has proved especially adept at handling the Burgers equation 
(Shu 1999). Recently, a hybrid hydrodynamic/A^-body code based on the WENO scheme 
has been developed. It has passed typical tests, including the Sedov blast wave and the for- 
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mation of the Zel'dovich pancake (Feng et al. 2004). This code has been successfully used to 
produce the quasi-stellar object (QSO) Lya transmitted flux, including the high-resolution 
sample HP1700+6416 (Feng et al. 2003). The statistical features of these samples are in 
good agreement with observed features not only on second-order measures, like the power 
spectrum, but also up to orders as high as 8 for the intermittent behavior. Hence, we believe 
that it would be reasonable to study the entropy and temperature fields of the IGM with 
the simulation samples given by the WENO scheme. 



3.2. Samples 

For the purpose of this paper, we run the hybrid A^-body/hydrodynamic code to trace 
the cosmic evolution of the coupled system of both dark matter and baryonic gas in a flat low- 
density ACDM model, which is specified by the cosmological parameters {flm, ^a, h, erg, fib) = 
(0.3,0.7,0.7,0.9,0.026). The baryon fraction is fixed using the constraint from primordial 
nucleosynthesis as Qj, — 0.0125/i~^ (Walker et al. 1991). The linear power spectrum is taken 
from the fitting formulae presented by Eisenstein & Hu (1998). 

The simulations were performed in a periodic, cubic box of size 25 h~^Mpc with a 192^ 
grid and an equal number of dark matter particles. The simulations start at a redshift 
z — 49, and the results are noted at the redshifts of z =4.0, 3.0, 2.0, 1.0, 0.5, and 0.0. The 
time step is chosen by the minimum value among the following three timescales. The first 
is from the Courant condition given by 

cfl X a{t)Ax 

< 71 \ i i \ \ r 

max{\v^\ + Cs, \Vy\ + Cg, \v^\ + c^) 

where Ax is the cell size, Cg is the local sound speed, Vx, Vy and Vz arc the local fluid velocities, 
and cfl is the Courant number, for which we take cfl = 0.6. The second timescale is imposed 
by cosmic expansion, which requires that Aa/a < 0.02 within a single time step. The last 
timescale comes from the requirement that a particle moves not more than a fixed fraction 
of the cell size. 

We then produced two realizations. One includes the radiative cooling and heating, and 
one does not. For the latter, the initial temperature of gas is taken to be 10^ K. For the for- 
mer, the atomic processes, including ionization, radiative cooling, and heating, are modeled 
in the same way as in Cen (1992) in a plasma of hydrogen and helium of primordial compo- 
sition [X = 0.76, Y — 0.24). The processes such as star formation, feedback due to SNs, and 
AGN activities have not yet been taken into account. A uniform UV background of ionizing 
photons is assumed to have a power-law spectrum of the form J(i/) = J21 x 10~^^ {i'/i/Hi)~"erg 



(11) 
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s^^cm^^sr^^Hz"^, where the photoionizing flux is normahzed by the parameter J21 at the 
Lyman hmit frequency vhi and is suddenly switched on at 2; > 10 to heat the gas and 
reionize the universe. 

For statistical studies, we randomly sampled 500 one-dimensional (ID) fields from the 
simulation results at every redshift of z=4, 3, 2, 1, 0.5, and 0.0. Each one-dimensional sample, 
of size L=25 h~^Mpc, contains 192 data points. For each of these points, information about 
the mass density and peculiar velocity of dark matter, and mass density, peculiar velocity, 
and temperature of the IGM are recorded. All of the statistical analyses below are based on 
the first realization, i.e., the one including cooling and heating. The second one is only used 
to produce Figure 9, which shows the effects of cooling and heating. 

3.3. An example of the IGM field 

We select one sample as an example of the IGM field and show in Figure 1 the one- 
dimensional spatial distributions of the temperature, mass density, and peculiar velocity of 
the IGM at the two redshifts z={) and 2;=0.5. For comparison, we also show the densities 
and velocities of dark matter in the corresponding panels. For identification, we use the 
notations pdm and pigm for the mass densities of the dark matter and IGM, respectively (p 
in §2 actually is Pigm-) The mass densities pdm and pigm are in units of and Pj^^, which 
are the mean densities of dark matter and gas, respectively. Spatial position is described by 
the comoving coordinate x. 

Figure 1 shows clearly that the relations between the density and velocity fields of the 
dark matter and IGM are very complex. The density and velocity fields of the IGM are not 
always faithful tracers of the underlying dark matter fields, although in the linear regime 
they are point-by-point tracers on scales larger than the Jeans length (Bi et al. 1992; Fang 
et al. 1993; Bi 1993; Nusser 2000; Matarrese & Mohayaee 2002). Many of the discrepancies 
between the IGM and dark matter fields shown in Figure 1 can be explained by the Jeans 
smoothing (Theuns et al. 2000). However, Figure 1 shows that several differences between 
the mass fields of the IGM and dark matter are on scales of a few h^^ Mpc, which is larger 
than the corresponding Jeans length of the IGM. 

The big jump in the IGM temperature around the position a; ^ 13 hr^ Mpc is typically 
caused by a strong shock. The temperature of the preshocked gas is ~10'^K and increases to 
lO^K by the shock. In this case, the temperature increases by a factor of ~ 10^. According 
to the shock theory of a polytropic gas (Landau & Lifshitz 1959), the Mach number should 
be M ~ VlO^ ~ 30. The value is reasonable. Moreover, for strong shocks, the density has 
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to increase by a factor of ~ (7 + l)/(7 — 1) =4. This is also consistent with the density 
field shown in Figure 1. It is also interesting to note that the jump exists in the region of 
Pdm < 1) i-G., not located close to filaments of dark matter. Therefore, gravitational shocks 
of the IGM can develop not only around massive dark halos but also in low-density regions. 

4. Temperature field of the IGM 

4.1. Multiple phases in the temperature-density diagram 

In Figure 2 we show the relations between the temperature and mass density of the IGM 
at z =4, 3, 2, 1, 0.5, and 0.0. Each panel in Figure 2 contains 19,200 data points randomly 
drawn from a total of 500 samples at each redshift. The basic feature of the T — pig^ 
diagram looks similar to those given by other cosmological hydrodynamic simulations or 
semi-analytical models (e.g., Ryu et al. 1993; Katz et al. 1996; Theuns et al. 1998; Dave et 
al. 1999; Valageas et al. 2002; Springel &: Hernquist 2002). Roughly speaking, there are two 
phases in the T — pig^ plane. (1) The first phase manifests a tight correlation in the bottom 
left regions T (< 10^'^K) and 0.01 < pigm < 2, and (2) the second phase lies above the tight 
correlation regions and has a scattered distribution in which the temperature generally is 
greater than lO^K. 

The tight correlation lines can be well approximated by a power law as T cx with 
a slightly increasing from ~0.61 at z=A to ~0.69 at z=0, which is the same as the so-called 
IGM equation of state given by Hui & Gnedin (1997). The small discrepancy of a from the 
adiabatic index 7 — 1 = 2/3 results from the photoionization heating and radiative cooling 
(Valageas et al. 2002). 

The gas in phase 2 sometimes is called the warm- hot intergalactic medium (WHIM), as 
its temperature is ~ 10"'"^ — 10'' K, i.e., higher than the gas in phase 1 by a factor of about 
10 - 10^ (Gen & Ostriker 1999). From Figure 2 we can see that the WHIM underwent a 
significant evolution from redshift z — Ato the present. At 2; = 4 WHIM events are rare and 
mostly reside in the regions of mass density, pig^ > 1, while WHIM becomes quite common 
when z <2. Note that the high-temperature (10^ — 10^ K) events, and hence strong shocks, 
can also occur in the regions with density pigm ~ 1, and even pig^ — 0.1. In other words, 
the case shown in Figure 1 is common. 

Comparing with numerical results cited above by other groups, phase 2 gas is more 
significant in our simulation. Especially in the low mass density {pigm < 1) regions, the 
phase 2 in Figure 2 is pronounced. This is because the WENO code is more effective in 
capturing shocks in both high- and low-density regions. 
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4.2. Statistical decoupling between the IGM and dark matter 

From Figure 1 wc have already seen the discrepancy between the IGM and dark matter 
fields. The Pigm/ Pdm ratios are not equal to unity everywhere. Hence, the baryonic fraction 
is not always equal to the cosmological average value = VLi^/VLrn = 0.085. To show this 
deviation, we plot the relationships between the mean baryonic fraction and temperature in 
Figure 3. It is obvious that the baryon fraction at 2; < 1 significantly exceeds the cosmological 
value fb for the WHIM (at lO^K < T < 10^ K) . This is a well-known result nowadays (Cen 
& Ostriker 1999; Dave et al. 2001). The concentration of baryons into WHIM can be seen as 
follows. The gas is compressed by a factor of ~ 4 by strong shocks. On the other hand, the 
dark matter mass density at the pre- and postshock is almost the same. Thus, the baryon 
fraction in the postshocked IGM can be enhanced as high as 4 x 0.085 ~ 0.3. This result is 
consistent with what is shown in Figure 3. 

However, WHIM may not necessarily reside in dense regions. In Figure 4 we present a 
cumulative probability distribution P{> Pdm) foi^ WHIM (T > 10^ K) at different redshifts, 
which is the fraction of WHIM with dark matter density larger than pdm- Note that we use 
the dark matter density pdm in Figure 4, not Pigm, as massive dark halos should be identified 
by high pdm, rather than by high pigm- If gravitational collapsing or collapsed regions are 
characterized by pdm > 10, Figure 4 shows that most of WHIM at 2; < 2 does not reside 
in or close to massive dark halos. In other words, many strong shocks may not be due to 
the accretion by the filaments of dark matter. One should not simply identify all WHIM as 
dense clouds or massive halos. 



4.3. Probability Distribution Fianctions of temperature field 

The multiphased features of the temperature field can be described by the PDF of the 
temperature field. The results are shown in Figure 5. We use three density ranges: (1) low, 
0.1 < Pdm < 0.4, (2) medium, 0.4 < pdm < 2.0, and (3) high, 2 < pdm < 10. All PDFs have 
similar shape but shift to higher temperature with higher density. The PDFs are sharply 
peaked at some low temperatures, while some have long tails on the high-temperature side. 
The tails are approximately given by a power law T~" with the index a ~ 4 at 2; = 4, and 
a ~ 1.7 at z = 0. Therefore, the high-temperature tails are prolonged at low redshifts. 
A field with such power-law-tailed PDF is typically intermittent (Nakao 1998). We also 
calculate the T — pdm diagram, which shows the similar result of multiphased structure as 
the T — Pigm diagram. 

We show in Figure 6 the relations between the mean temperature T and the dark matter 
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mass density pdm- One can see that T is actually inadequate to describe the temperature 
field if the field is long-tailed, or multiphascd. From Figure 6 we can see that T at redshifts 
< 1 and ~ 1 is in the range 10'^. 5 — 10^ K. However, Figure 5 shows that in the range 
'Pdm ^ ^ probability for temperature of 10'^ — lO'^K is comparable to the probability for 
T = 10^ — 10^ K. Therefore, the IGM cannot be approximately modeled as a single phase 
gas with temperature of either ~ 10^.5 — 10^ K or 10'^ — lO'^K. 

In Figure 6 each panel has three T — pdm curves, which correspond to the data smoothed 
by windows on scales 0.26, 1.04 and 4.17 h^^ Mpc, respectively It is interesting to see that 
the T — Pdm relations arc basically independent of the smoothing scales. Therefore, the 
multiphascd feature of the IGM temperature does not come from the sampling, but from 
the non-Gaussian PDF of the field. Such a field is substantially different from a uniform or 
a Gaussian field. 

The T — pd^ relation of WHIM given by Dave et al. (2001) is approximately T oc Pigm- 
From Figure 6 we can also see T oc pig^- However, they are different. The mean T in Figure 
6 is given by averaging over all pixels including both the WHIM and the gas with low (< 10^ 
K) temperature, since the T — pdm relation for low (< 10^ K) temperatures is T oc plg^ with 
a < 1. Therefore, for WHIM only, we may have T oc p"^^ and a > 1. Therefore, Figure 6 
shows again that our samples contain more shock-heated events. 



5. Entropy field of the IGM 

5.1. Entropy-temperature and entropy-mass density relations 

The entropy of the IGM is measured by the parameter T/n^^^^, where Uigm is the number 
density of baryons. It is actually related to the entropy density s of the IGM by s/ pigm = 
\nT / nfg^+ const. As we will not use s below, we use this notation for the entropy parameter 
or simply entropy, S = T/n^g'^. If temperature T is in units of keV (1 keV ~ 1.2 x 10^ K), 
entropy S is in units of keV cm^. Figure 7 presents the relation between the entropy and 
temperature of the IGM at z = 4, 3, 2, 1, 0.5, and for the 19,200 data randomly drawn 
from the total 192^ data points. 

Figure 7 also shows the two phases of the IGM: (1) a tight correlation along a line 
of ~ const at the bottom of the figure, i.e., the entropy of the IGM is independent of 
IGM temperature, and (2) a scattered distribution above the tight correlation fine. They 
are the same as Figure 2. The gas in phase 1 approximately satisfies T oc pfg^- The gas in 
phase 2 experiences an entropy increase process, mainly because of the heating by shocks 
and complex structures in the IGM. A strong shock increases T by a factor of 10^ — 10^, 
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iT-igm by a factor of ~4, and therefore, the postshocked entropy can be enhanced by a factor 
of 10^ - 10^ 

Phase 2 shows a significant evolution from redshift z — A to the present. There are 
almost no high-entropy {S > 50 h~^/^ keV cm^) events when z — A. At z — 2 and 3, high- 
entropy events are biased to T > 5 x 10~^ keV regions. When z < 1, high-entropy events 
occur in all temperatures T > 10~^ keV, or greater than 10^ K ranges. 

The S - pdm diagram is shown in Figure 8. We again use pdmi not pigmi as the clustering 
of the mass field is indicated by high-p,^^. Although pj^^ does not trace p^m point by point, 
one can still see the two phases: a dense distribution along a line of constant and a 
scattered distribution above the line, given by the entropy increase because of shocks. For 
redshift 2; < 3, the gas in phase 2 can reside in both low {pdm < 1) and high {pdm > 1) dark 
matter density areas. 

Figure 9a presents the mean entropy S versus density pdm- The entropy of the IGM 
increases by about 2 orders of magnitude from redshift z = 4 to 0.5. In the redshift range 
3 > 2; > 1, the mean entropy S is higher in higher density areas. For dark matter mass 
density pdm > 10^, the mean entropy S is greater than 60 h~^/^ keV cm^ at the epoch z ~ 2, 
and greater than 20 h~^/^ keV cm^ at ;2 ~ 3. At 2; < 0.5, the shocks and complicated 
structures yield an "entropy floor" S ~ 100 — 300 h~^^^ keV cm^ a,t z < 0.5 over the entire 
IGM from low-density {pdm — 10~^) areas to high-density {pdm — 10^) regions. This is the 
cosmological background of entropy floor generated by the nonlinear evolution of the IGM 
in the gravity of dark matter. 

In Figure 9b we show the S - pigm relation for samples produced by the same cosmological 
parameters, but without considering radiative cooling and photoionization heating. We see 
that the relation S - pigm at 2; < 2 in Figure 9b is about the same as that of Figure 9a. 
Especially in the high-density range {pdm > 10), Figures 9a and 9b are the same. Therefore, 
the entropy increase of the IGM at z < 2 shown in Figure 9a is not due to the heating 
and cooling, but mainly determined by shocks and complicated structures. The effects of 
radiative cooling and heating to the entropy evolution are negligible at z <2 and high-density 
areas. At higher redshift z > 2 and low-density areas, however, the effects of radiative cooling 
and heating become important. 

It is interesting to compare our result of the entropy floor with the semi-analytical 
estimation by Valageas et al. (2003, hereafter VSS03). Using the data of Dave et al (2001), 
they found that WHIM can provide an entropy floor (or the mean S) of the order of 10^ h~^^^ 
keV cm^ at 2; < 0.5, which is about the same as Figure 9. However, their result is given by 
an average of S over all WHIM. As has been pointed out in §3.2, most WHIM at 2; < 2 does 
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not reside in high-density areas. Therefore, it is unclear whether the high-entropy floor given 
by an overall average of the WHIM is in the high-density areas. Moreover, the average of S 
in VSS03 is done by the method of adiabatically bringing all patches into a same effective 
density. With this method, the low-density areas would have significant contribution to the 
mean S. In our study, the mean of S is taken in each given bin of the dark matter density, 
and therefore from Figure 9 one can unambiguously conclude that the high-entropy floor of 
~ 10^ h~^/^ keV cm^ first formed in high-density p^m > 10 areas since z — 2, i.e., the areas 
that are favored by the formation of galaxy clusters. 

Moreover, the mean of S given by VSS03 falls steeply with redshift. VSS03 found 
5* < 10 h^^^'^ keV cm^ at z > 1.5. Figure 9 also shows the decrease of S with redshift. 
However, it is much slower than that of the VSS03. Figure 9 shows that the entropy fioor is 
still ^ > 10 h~^/^ keV cm^ when z — 3 and pdm > 10. The difference between the redshift 
dependence of the VSS03 and Figure 9 is caused by two reasons. (1) The number of S of the 
VSS03 is given by an average over all density regions, and Figure 9 shows that the decrease 
of S with redshift is faster in lower density {pdm < 1) regions; and (2) the WENO code is 
more effective in describing shocks than the data used in VSS03. 

5.2. PDF of entropy field 

At first glance, the mean entropy of Figure 9 may contradict the Lya absorption in 
QSO's absorption spectrum. A uniform entropy floor S ~ 100 h~^/^ keV cm^ corresponds 
to T ~ 1.8 X 10^ K at 2; = 1 and pigm — 1- This temperature is obviously unacceptable 
with the observations of Lya forest lines. Even for a uniform entropy floor as small as 10 
h~^^^ kcV cm^, the IGM temperature is about 4x10^ K at 2; = 2 and pig^ — 1- This is 
also contradictory with the small-scale intermittent spikes in the transmitted flux of QSO's 
LycK absorption. This argument is often used to constrain preheating models. However, this 
argument no longer holds if the IGM is multiphased. This can be seen from the PDF of the 
entropy field. 

Similar to the temperature field, the entropy field is also highly non-Gaussian. Figure 
10 presents the PDFs of the entropy field at redshifts from 2; = 4 to 0. The shapes of the 
PDFs at different redshifts are similar, but the curves shift to higher S at lower z. The 
PDFs are peaked at the low-entropy end and gradually decrease with the increasing 5" by a 
power law oc S~'^ with 6 ~ 4 at 2; ~ 4 and b ~ 1.6 at 2; ~ 0. That is, the tail is longer when 
the redshift z decreases. Therefore, in terms of entropy, the IGM is also in multiple phases. 
High- and low-entropy gas coexist in the IGM. In this case, the mean entropy, S, does not 
mean uniform heating. For a long-tailed PDF, the mean is much larger than the median. 
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Thus, most of the space actually has entropy S less than the mean. 

Figure 11 gives the volume filling factor V{> S), which is the volume fraction with 
entropy larger than a given S. It shows that the volume fraction of high entropy {S > 
100/i~^/^ keV cm^) is very small for all redshifts. Although the mean entropy S is as high as 
~ 200/i-^/=* keV cm^ at ^ = (Figure 9a), the volume fraction V{> 200) is only about 10%. 
Therefore, it is not surprising that the high- and low-entropy phases can coexist in the IGM. 
To approximate such a multiphased field by a uniform distribution would be misleading. 

5.3. Mass fraction of high-entropy and high-density regions 

We now study the distribution of entropy of the IGM in mass space. We are interested in 
the regions with high entropy and high density, as they provide the environment of baryonic 
structure formation. Let us consider the regions with high dark matter density p^m = 5 or 
10 to be the hosts of collapsed baryonic clumps. Figure 12 shows the baryonic mass fraction 
M(> S) of regions with the entropy larger than a given S and dark matter density pdm > 5 
and > 10, respectively. 

First, we see from Figure 12 that M(> 10) ~ M(> 100) at ^ = for either pdm > 5 or 
> 10. That is, at z = 0, no high-density regions pdm > 5 have entropy less than 100 h~^^^ 
keV cm^. At higher z, the mass fraction with higher entropy becomes smaller. Nevertheless, 
there are already ~ 10% baryons at 2; ~ 2 to have entropy greater than 50 h~^^^ keV cm^. 
This feature is not sensitive to whether pdm > 5 or 10. Since the baryon fraction in collapsed 
halos is about 10%, we may conclude that all structure formation is probably in the entropy 
environment of 5" > 50 h~^^^ keV cm^ if the gas falls into the halos not earlier than 2; ~ 2. 

In calculating Figure 12, the mass field pdm was decomposed to the size of the resolution 
of the simulation, which is 0.13 Mpc. To study the scale dependence of the entropy 
floor, we decomposed the dark matter mass held and entropy field into different scales. This 
decomposition can be done by the orthogonal scaling functions of discrete wavelet transform 
(DWT). As has been showed (Xu et al. 1998), the high-density regions on the scale 1.5 h"^ 
Mpc given by the DWT decomposition can be identified as clusters, which yields the same 
result as traditional identification schemes, such as the friends-of-friends algorithm. 

Figure 13 presents the relations between the mean entropy S and pdm for the mass field 
decomposed by the DWT scaling function on scales L = 0.26, 1.04 and 4.17 Mpc. The 
mean S in Figure 13 is given by averaging over all regions that are on the given scale and 
have the dark matter density pdm in a given bin. Therefore, S* at L = 1.04 h^^ and pdm > 10 
is the mean entropy of the IGM in the regions that host the formation of groups and clusters. 
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Figure 13 also shows that the S-pdm relations are basically scale- independent when 
z >?>. At small rcdshifts -2 < 2, however, one can see that the larger the size L, the higher 
the mean entropy is in high-density regions. That is, the high-entropy IGM is biased to the 
places that have high dark matter density in a large size area. This is probably because 
more shocks occur around massive halos. Thus, at epochs -2 < 2, massive halos are generally 
surrounded by high-entropy baryonic gas. This result shows again that the entropy floor in 
the regions with high mass density and on large scales is generally larger than the estimation 
by an over-all average (VSS03). 



6. Discussions and Conclusions 

We have studied the evolution of the temperature and entropy fields of the IGM in 
the nonlinear regime of gravitational clustering. We have used the iV-body /hydro dynamic 
simulation samples produced by the WENO code, which effectively describes shocks. Strong 
shocks can significantly change the mass density and velocity of the IGM fiuid, at the same 
time having no effect on the dark matter field. Therefore, when the shocks emerge, the 
dynamical similarity between the IGM and dark matter is broken. The gravitational shocks 
and complex structures of the IGM can take place in nearby massive dark halos as well as 
in low mass density regions. 

The entropy production by shocks makes the IGM multiphased. The temperature T and 
entropy S of the IGM arc no longer related to the mass density pigm or pdm by a single- value 
relation. The IGM contains gas in both high- and low-entropy phases, or both high- and 
low-temperature phases. For a given pig^ or p^^, the IGM temperature and entropy lie in 
very wide ranges from ~10'^ to ~10^ K, and from ~1 to ~10'^ h^^^^ keV cm^, respectively. 
Therefore, the temperature and entropy fields of the IGM cannot be approximately described 
by a randomly uniform Gaussian distribution. 

We should point out that star formation and their feedback on the IGM evolution 
are not considered in our simulation. Roughly, there arc two types of the feedbacks: (1) 
photoionization heating by the UV emission of stars and AGNs, and (2) injection of hot 
gas and energy by SN explosions or other sources of cosmic rays. The photoionization 
heating actually can be properly considered if the UV background is adjusted by fitting the 
simulation with observed mean fiux decrement of QSO's Lyct absorption spectrum (Feng et 
al. 2003). The effect of injecting hot gas and energy by SNs is localized in massive halos. 
For instance, metal absorption systems (like CIV and Mgll) in QSO spectra are generally 
associated with collapsed halos (e.g., Bi & Fang 1996). A recent Lya observation around 
protoclusters (Adelberger et al. 2003) implies that AGN heating does not drastically affect 
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the gas in clusters. Nor does it affect tfie IGM in low-density areas. Therefore, these heating 
mechanisms are probably not strong enough to change the basic feature of the multiphased 
temperature and entropy field given in this paper. 

The multiphased feature and non-Gaussianity of the IGM temperature and entropy 
fields can effectively be apphed to explain the following IGM-related observations. 

1. Entropy floor of groups and clusters. — For all regions on scales greater than 1 
Mpc and with dark matter mass density pd^n > 10^, the mean entropy S, or entropy floor, 
is greater than 80 h"^/^ keV cm^ since the epoch 2; ~ 2 and greater than 20 h~^^^ keV cm^ 
at 2; ~ 3. Even for regions of pdm > 10, the mean entropy still have ^ > 25 h~^^^ keV 
cm^ at the epoch 2; ~ 2 and 10 h~^^^ keV cm^ at 2; ~ 3 (Fig. 13). Therefore, gravitational 
shock is an important provider of high-entropy gas around massive halos, which are the 
hosts of clusters and groups. This result is adequate to explain the entropy excess from the 
observations of clusters and groups if the epoch of the gas falling in cluster cores is not earlier 
than 2; ~ 2 — 3 (Ponman et al. 2003). In other words, the IGM is already significantly 
preheated by gravitational shocks at the epoch 2; ~ 2 — 3, and therefore other sources of 
heating may not be important if the cluster core with entropy excess formed after that epoch 
(Kay & Bower 1999; Borgani et al 2002). 

2. LycK forest hues. — The Lya forest lines in QSO's absorption spectrum of 2 < 2; < 4 
have a column density Nhi of HI atoms in the range 13 < In Nhi < 14, which corresponds 
to areas of pdm — 1 (Bi & Davidsen 1997). The thermal broadening of these hues shows 
that the temperature of the IGM at this density area is about 10^ - lO^K. The high-entropy 

floor (S ~ 100h~^^^ keV cm^) will contradict the Lya forest if the entropy field is uniform. 
However, the temperature and entropy fields are highly non-Gaussian and of multiple phases. 
The volume of the universe is actually dominated by regions of low temperature and low 
mass density. For all redshifts 2; < 4, more than 90% volume with dark matter mass density 
{Pdm < 2) is occupied by the IGM with temperature less than lO'^'^K (Fig. 5). This provides 
the room for the QSO's Lya forests. The mean temperature is higher for lower redshifts 
(Fig. 6). Therefore, the number density of the forests lines, or absorbers, is lower at lower 
redshifts. These properties match the constraints on the IGM from QSO's Lya forests. 
Moreover, the Lya forests are very well modeled by the semianalytical lognormal model, 
which assumes that the PDF of the IGM mass density field is lognormal, and long-tailed (Bi 
& Davidsen 1997). This assumption basically is consistent with the long-tailed PDF of IGM 
temperature field (Fig. 5). 

3. Intermittency of Lya transmitted flux. — The Lya transmitted flux of the QSO's 
absorption spectrum is found to be remarkably intermittent (Jamkhedkar et al. 2002, 2003). 
That is, the power of the flux fluctuations concentrates in the rare spikes. The intermittent 
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features can be detected even on scales corresponding to thermal broadening of gas with 
T ~ 10^ — 10^ K. Different from Lya forests, the intermittency of the Lya transmitted 
flux does not become weaker, but become even stronger at lower redshifts. It is impossible 
to explain the intermittency with single-phase IGM with temperature ~ 10^ K. However, 
strong shocks yield many big discontinuities of the IGM temperature (Fig. 1), and therefore 
the power of the temperature fluctuations is also spiky. The long-tailed PDFs also reflect 
the intermittency of the temperature field. The lower the redshifts, the longer the tails, 
and therefore the intermittent features are more prominent at lower redshifts. Moreover, 
the IGM at density pdm — 1 regions contains components of T ~ 10^ — 10^ K as well as of 
T ~ 10^ — 10^ K. Therefore, the intermittency of the Lya transmitted fiux can be well fitted 
with these simulation samples (Pando et al. 2002; Feng et al. 2003). 

All these results are in good agreement with the theoretical expectation of §2. In the 
nonlinear regime, the Reynolds number is larger and leads to the Burgers turbulence. There- 
fore, the gravitational shock in the IGM provides a unified explanation of the features of the 
IGM in both low- and higher density regions. It includes the departure of the IGM-dark mat- 
ter similarity mapping, formation of the entropy excess groups and clusters, the temperature 
and entropy in the regions of Lya forests and the intermittency of Lya transmitted flux. 
The IGM-dark matter discrepancy gives a deep understanding of the nonlinear evolution of 
the system of the IGM plus dark matter. 
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Fig. 1. — One-dimensional spatial distributions of temperature T, mass density Pigm, and 
velocity of the IGM (solid lines), and mass density pdm and Vdm of dark matter (dashed lines) 
at z—0 and 0.5. The densities pdm and pigm are in units of and p^g^, which are the mean 
densities of dark matter and gas, respectively. 

Fig. 2. — Temperature T vs. IGM density pig^. Each panel is given by 19,200 points 
randomly drawn from the simulation sample in 25^h~^ Mpc^ box with 192^ data points. 

Fig. 3. — Relations of the mean baryon fraction with respect to temperature of the IGM. 
The redshifts are taken to be 4.0, 2.0, 1.0, and 0. 

Fig. 4. — Relations of the cumulated probability of WHIM (T > 10^ K) with respect to the 
dark matter density pdm- The redshifts are taken to be 4.0, 2.0, 1.0, and 0. 

Fig. 5. — PDFs of the IGM temperature field in the ranges of 0.1 < pdm < 0.4, 0.4 < pdm < 2, 
and 2 < pdm < 10, respectively. The redshifts are taken to be 4.0, 2.0, 1.0, and 0. 

Fig. 6. — Mean temperature T vs. pdm for samples smoothed on scales Lj = 0.26, 1.04, and 
4.17 h"-*^ Mpc, respectively. The redshifts are taken to be 4.0, 2.0, 1.0, and 0. 

Fig. 7. — Entropy S vs. Temperature T. Each panel is given by 19,200 points randomly 
drawn from the simulation sample in a 25^/i-3 Mpc^ box with 192^ data points. 

Fig. 8. — Entropy S vs. dark matter density pdm- Each panel is given by 19,200 points 
randomly drawn from the simulation sample in a 25^h~^ Mpc^ box with 192^ data points. 

Fig. 9. — Mean entropy S vs. pdm (a) for the simulation with ionizing, radiation heating, 
and cooling at redshifts 4.0, 3.0, 2.0, 1.0, 0.5, and 0.0, and (b) for the simulation without 
ionizing, radiation heating, and cooling at redshifts 4.0, 2.0, 0.5, and 0.0. 

Fig. 10.— PDFs of entropy field at redshifts 4.0, 3.0, 2.0, 1.0, 0.5, and 0.0, respectively. 

Fig. 11. — Volume filling factors of entropy S field at redshifts 4.0, 3.0, 2.0, 1.0, 0.5, and 0.0, 
respectively. 

Fig. 12. — Baryonic mass fraction with entropy larger than 10, 50, 100, 200, 500, and 1000 
/i^^/^ keV cm^, and pdm larger than 5(lcft panel) and 10 (right panel), respectively. 

Fig. 13. — Mean entropy S vs. pdm for samples smoothed on scales Lj — 0.26, 1.04, and 
4.17 h-i Mpc. The redshifts are taken to be 4.0, 3.0, 2.0, 1.0, 0.5, and 0. 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405 1 39v2 




T(K) 



10" 



10"^ 



10" 



Oh 10 



10" 



10" 



10" 



I I 1 1 1 lll 1 — I I I 1 1 llj 1 — I I I I lllj 1 — I I I 1 1 lll 1 — I I I 1 1 

z=4 

0.1<Pdm<0.4 

0.4<pdm<2.0 
2.0<pd^<10.0 




I 1 1 iiii i f 



z=l 

0.1<Pdm<0.4 

0.4<pdm<2.0 
2.0<pdm<10.0 



jj I 



I 1 1 1 ii| 1 — I I 1 1 1 llj 1 — I I I I lllj 1 — I I I 1 1 ii| 1 — I I I 1 1 ii| 1 — I I 1 1 1 II 

z=2 

0.1<Pdm<0.4 

0.4<pdm<2.0 
2.0<pdm<10.0 




z=0 

0.1<Pdm<0.4 
0.4<Pdm<2.0 

2.0<pdm<10.0 




I I I M il l 



10^ 10^ lO"^ 10^ 10^ 10^ 10^ 10^ lO"^ 10^ 10^ 10^ 10^ 



T(K) 



T(K) 



This figure "fig7.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405 1 39v2 



This figure "fig8.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405 1 39v2 



This figure "figl3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405 1 39v2 



